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ABSTRACT 

Using light curves and host galaxy spectra of 101 Type la supernovae (SNe la) with redshift z < 0.3 
from the SDSS Supernova Survey (SDSS-SN), we derive the SN la rate as a function of progenitor 
age (the delay time distribution, or DTD). We use the VESPA stellar population synthesis algorithm 
to analyze the SDSS spectra of all galaxies in the field searched by SDSS-SN, giving us a reference 
sample of 77,000 galaxies for our SN la hosts. Our method does not assume any a priori shape for the 
DTD and therefore is minimally parametric. We present the DTD in physical units for high stretch 
(luminous, slow declining) and low stretch (subluminous, fast declining) supernovae in three progenitor 
age bins. We find strong evidence of two progenitor channels: one that produces high stretch SNe la 
< 400 Myr after the birth of the progenitor system, and one that produces low stretch SNe la with 
a delay > 2.4 Gyr. We find that each channel contributes roughly half of the Type la rate in our 
reference sample. We also construct the average spectra of high stretch and low stretch SN la host 
galaxies, and find that the difference of these spectra looks like a main sequence B star with nebular 
emission lines indicative of star formation. This supports our finding that there are two populations 
of SNe la, and indicates that the progenitors of high stretch SNe are at the least associated with very 
recent star formation in the last few tens of Myr. Our results provide valuable constraints for models 
of Type la progenitors and may help improve the calibration of SNe la as standard candles. 
Subject headings: Cosmology: Distance Scale, Galaxies: Stellar Content, Stars: Supernovae: General 



1. INTRODUCTION 

Type la supernovae (SNe la) are thought to be the 
products of a thermonuclear explosion of a white dwarf 
reaching the Chandrasekhar mass - the mass limit al- 
lowed by the supporting electron degeneracy pressure. 
This requires the white dwarf to accrete mass from a 
binary companion, though the nature of this companion 
and the timescale of the accretion are hotly debated. The 
single degenerate model proposes the companion to be a 
main sequence or giant star with mass transfer by Roche 
lobe overfiow, while the double degenerate scenario pro- 
poses that the companion is also a white dwarf and the 
two objects merge. While the connection between these 
two scenarios and the delay time from progenitor birth 
to SN la explosion is unclear, knowledge of this delay 
would constrain the initial masses of SN la progenitors. 

Interest in SNe la has greatly increased since their 
use as standardized candles led the di scovery of the ac- 
celerated expansion of the Universe (|Riess et al.l 119981 : 
iPerlmutter et al.lll999( ). The latest generat ion of SN la 
surveys, includi ng SNLS (|Astier et al. 2006), ESSENCE 
(IWood-Vasev et al. 20071) . and SDSS-SN (Kessler ct all 
I2009D . however, are limited by the possibility of system- 
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atic uncertainties in their calibration of SN la luminosi- 
ties. This has sparked great interest in any physical prop- 
erties that coul d alter the p eak luminosity/light curve 
width relation (|Phillipslll993D used to calilDrate SNe la. 
In particular, it has highlighted our incomplete knowl- 
edge of SN la progenitor systems. 

Due to their low luminosities, SN la progenitors 
have never been directly observed. Constraints on 
their nature and rates of explosion must therefore come 
from studies of SN la environments, and a variety of 
techniques ha.ye been used to carry out s uch studies. 
iSullivan et"all (|2006D . fMannucci etldl (I2006D. and others 
use host galaxy photometry and [Gallagher et al.l (|2005[ ) 
use spectra to derive stellar populations and metallici- 
ties, and correlate these luminosity- weight ed average val- 
ues ag ainst SN la ra tes and properties. ISullivan et al.l 
and Ma nnucci et al.l s ee higher SN la rates in star- 
forming galaxies, while iGallagher et al.l find that spirals 
host more luminous SNe. To de rive constraints o n the 
ages of the proge n itor s ystems, INeill et all (j2006[) and 
IGal-Yam fc Maol (I200I compare the evolution of the 
supernova rate with the evolution of the cosmic star 
formation rate. Because this rate changes slowly with 
redshift in the local universe, the metho d places only 
weak constraints on the progenitor ages. ICooper et al.l 
(|2009D have used host galaxy clustering as a proxy for 
the metallicities of SN la environments and found a sig- 
nificant correlation with supernova rate or luminosity. 
Badenes et al. (2009) have measured the /omZ stellar pop- 
ulations for four historical SNe la in the Large Magellanic 
Cloud, finding that three of the supernovae live in regions 
with little star formation over the last several Gyr. 

Over the past decade and a half, analyses of SN 
la environments (including several papers cited above) 
have produced evidence that there are at least two dis- 
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tinct populations of SNe la: blue, star-forming galax- 
ies have higher supernova rates and host more luminous, 
slower-declining supernovae than do red, passive galaxies 
jHamuv ct al. 1996; Howell 2001; van den Bcrgh et alJ 
[200a .Mannuc ci et all [2005, but see ISchawin ski 2009). 
More recently, this conclusio n has been c o nfirm ed in 
large, local SN la samples by iHicken et al.l (|2009[ ) and 
Maoz et al.l (I2010D a nd in a higher redshift sample by 
Sullivan et al.l ()2010D . The observation of at least two 
populatio ns has le d to the so-called A + B model 
(jScannapi eco fc Bild stcn 2005), in which the supernova 
rate is modeled as the sum of a term proportional to total 
stellar mass (the delayed component) and a term propor- 
tional to recent star formation (the prompt component): 



SNR = AM^ + BAh 



(1) 



Values of A and B have been determined by 
ISullivan et all ()2006[ ). lNeiU et al.l (poM ) . iMannucci et al.l 
(1200511 and others, with significant scatter between the 
groups due to different SN la samples, methods for de- 
riving A and B, definitions o f stellar masses and p roxies 
for the star formation rates. lAubourg et al.l (|2008f ) used 
host galaxy spectra and found that the supernova rate 
per unit mass for young stars is ~ 500 times higher than 
for old stellar populations, a factor of ^ 5 higher than 
earlier results. They also found evidence that the time 
scale of the prompt component is < 180 Myr. 

Unfortunately, the A + B model provides only weak 
constraints on the ages associated with the prompt and 
delayed components. Because of the strong dependence 
of stellar evolution timescales on initial mass, more pre- 
cise ages can powerfully constrain the main-sequence 
masses of the progenitors. This will require a better ap- 
proximation to the full delay-time distribution (DTD), 
the explosion rate as a function of progenitor age. The 
DTD, e{t), has units of number of supernovae per unit 
stellar mass per year, and relates a galaxy's supernova 
rate to its star formation rate ip{t): 



SNR{t) ^ / e{t~t')i:{t')dt'. 



(2) 



Thus, e{t) represents the probability per year that a 
unit of stellar mass will produce a Type la supernova 
a time t after its formation. The DTD sets the rates and 
timescales of SN la production that must be matched by 
progenitor models. 

There have been several attemp t s to p arametrize and 
measure the DTD. iPritchet et all (|2008l) assume a con- 
tinuous, power-law DTD (e(t) cx t~P) and find that the 
SNLS survey constrains p to lie in the range 0.3 < p < 
0.7. Roughly speaking, this is required f or consistency 
with m easurements of A and B from, e.g.. ISullivan et aTl 
()2006( ). Using a theoretical argument to calculate the 
rate of white dwarf formation per unit stellar mass, they 
find that a single-degenerate model can yield a power 
law DTD of this form only if the fraction of white dwarfs 
which explode as SNe la is independent of progenitor 
mass. Because less m assive progenit ors should produce 
less massive remnants. IPritchet et al.l conclude that there 
must be another, pos s ibly d ouble degenerate, route to 
SNe la. iTotani et all (|2008D constrain the DTD in old 
stellar populations by selecting SN candidates in pas- 
sive galaxies from the Subaru/XMM- Newton Deep Sur- 



vey (SXDS, iFurusawa et~all I2008D . Combining these 
high redshift observations with the local observed SN 
la rate in elliptic als, they find e{t) cx in the range 
0.1 < t < 10 Gyr. ITotani et ahl argue that a DTD of this 
form supports a double-degenerate origin for delayed SNe 
la. 

In this paper, we aim to provide better constraints on 
the DTD and to determine whether luminous SNe la 
(with temporally "stretched" light curves) and less lumi- 
nous events have different progenitors from one another. 
Any attempt to measure the DTD relies on the intrinsic 
range of galaxy properties, and in particular, the differ- 
ences between supernova hosts and non-hosts. It is there- 
fore essential to have a supernova sample with a well- 
understood selection function and a well-defined control 
group of galaxies monitored for SNe. Otherwise, differ- 
ences between hosts and non-hosts may be attributed to 
biases or survey systematics rather than to supernova 
rates. Nearby surveys often comprise SNe discovered by 
many different techniques, making it difficult to define or 
build a control sample. High redshift searches are often 
blind (i.e. sensitive to SNe in all galaxies in the field), 
but lack spectra of the field galaxies. 

The Sloa n Digital Sky Survey Supernova S urvey 
(SDSS-SN, iFrieman et al]l2008l ISako et al.ll2008h sat- 
isfies both criteria. It is the only controlled, blind, dif- 
ference imaging supernova search at low redshift, with 
most supernovae in the range 0.05 ^ z < 0.35. SDSS 
(jYork et al.l 120001 ) also has a large spectroscopic cata- 
log of galaxies, of which about 83,000 lie in the stripe 
of sky monitored for supernovae. SDSS-SN is therefore 
the ideal survey with which to perform a statistical com- 
parison of the spectroscopic properties of hosts and non- 
hosts. 

For our study we select 101 SNe la from the SDSS- 
SN sample with acceptable light curves and SDSS host 
:alaxy spectra. We then use VErsatile SPectral Analysis 
Toieiro et"aLll2007D - VESPA - to derive star formation 
histories for all galaxies in the control sample. Given 
any assumed delay time distribution, these star forma- 
tion histories allow us to simulate a population of hosts. 
We then use the spectra themselves to compare these 
mock hosts to the observed SN la host galaxies. We 
compute likelihoods for a large number of assumed delay 
time distributions, thereby constraining the DTD. 

We have organized this paper as follows: in Section [2] 
we describe our datasets; in Sections |3| and ID we describe 
our methodology; in Section|5|we present our results and 
we finally discuss and conclude in Section [HI 

2. DATA: SDSS-SN 

The SDSS Supernova Survey (SDSS-SN ) was carried 
out along "Stripe 82" (jYork et al.1 l2000t ). a region of 
sky 2.5 degrees wide centered on the Celestial Equator, 
stretching from 20.7 hours to 3.9 hours (-50° - +59°) 
Right Ascension for a total area of 270 deg^. It was ob- 
served up to 80 t imes with the SDSS imaging camera 
(jCunn et al.lll998| ) ; supernova candidates were identified 
by difference imaging, and followed up with spectroscopy 
on oth er telescopes, and photometry from the SDSS tele- 
scope ()Gunn et al.ll2006D and othe r telescopes as well. 

S DSS-SN is describ ed in detail in lFrieman et al.l (|2008D 
and lSako et al.l ()2008[ ). Operating for three months of the 
year from 2005 to 2007, SDSS-SN has identified over 400 
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spectroscopically confirmed SNe la in the redshift range 
0.05 < z < 0.35, of which 146 fr om the 2005 seas o n hav e 
ugriz light curves pubhshed in iHoltzman et aLl ()2008D . 
Prehminary data for the remainder were kindly provided 
to us by the SDSS supernova team. We have selected 
definitive and probable SNe la this catalog of preliminary 
photometry using the types reported through the Central 
Bureau for Astronomical Telegrams'". We have further 
excluded two objects listed as peculiar SNe la, 2005hk 
and 2007qd. 

The preliminary photometric data lack errorbars, 
which we have estimated using the observed correlation 
between errors an d magnitudes in each b and in the pub- 
lished 2005 data ()Holtzman et al.l [20081) . Because our 
final analysis only uses the light curves to divide our 
sample into two broad subsamples, our results are insen- 
sitive to the details of these error estimates. We have also 
checked the reliability of the preliminary flux measure- 
ments using objects from the 2005 season for which final 
data are also available. The differences between prelimi- 
nary and final published fiuxes are approximately Cauchy 
distributed, with a mean offset from zero of ~ 0.5 /iJy 
and a full width at half maximum of ^ 2.5 /iJy. These 
are much fainter than a typical SDSS supernova peak 
magnitude of < 20, or a peak fiux of > 40 /xJy. 

The SDSS carried out spectroscopy of galaxies on 
Stripe 82, using the standard selection algorithms used 
throughout the SDSS survey. These include the Main 
Galaxy selection described in Strauss et al.' (2002), con- 
sisting of all galaxies with .Petrosian (1976) r-band mag- 
nitudes brighter than 17.77 and with a median redshift of 
or der 0.13, and the Lumin ous Red Galaxy (LRG) sample 
of lEisenstein et all (|200lD . which selects the most lumi- 
nous red ellipticals to z ^ 0.55. Because of the limited 
imaging carried out by SDSS in the Fall months, however, 
additional spectroscopy was carried out on Stripe 82 us- 
ing a variety of algori thms that went beyond the mai n 
spectroscopic sample (jAdelman-McCarthv et al.ll2006[ ). 
including a sample designed to calibrate photometric red- 
shifts, a sample of low-luminosity galaxies, an extension 
of the LRG sample to fainter and bluer objects, and a 
sample flux-limited in the u-band. In all, there are about 
83,000 unique galaxies on Stripe 82 with spectroscopic 
data, of which only 23,000 were selected using the main 
survey algorithms. We have further excluded all galax- 
ies with bad redshifts and, because we are interested in 
differences only in the stellar populations between SN 
la hosts and non-hosts, galaxies flagged as AGNs or 
QSOs by the SDSS pipeline. This leaves about 77,000 
unique galaxies from t he Seventh Data Release of SDSS 
(jAbazaiia n et al.ll2009D . of which 22,000 are in the Main 
Galaxy Sample. 

We have matched the SDSS-SN sample with this spec- 
troscopic sample, yielding a total of 133 SNe with host 
spectra in SDSS. The light curves of these supernovae 
range in quality, but 101 are sufficiently good to mea- 
sure parameters such as the peak flux and decline rate 
(see §3.21 and Figure [T] for examples). The 77,000 unique 
galaxies in Stripe 82 were not uniformly selected, and do 
not represent a flux-limited sample. However, our sub- 
sample of 101 SN hosts was selected from these 77,000 
galaxies based solely on the occurrence of a detectable 

^ http://www.cfa.harvard.edu/iau/lists/Supcrnovac.html 



SN la. These 101 SNe with SDSS hosts form the sample 
used in the rest of our analysis. 

3. METHODOLOGY: INITIAL DATA ANALYSIS 

To measure the delay time distribution, we need to be 
able to derive star formation histories for galaxies from 
their SDSS spectra. Because of the strong association 
of luminous (high stretch) SNe la w ith recent star for- 
mation seen by Sul livan et al.l (|2OO60 and others, we also 
investigate the dependence of the delay time distribution 
on supernova stretch. We begin by describing how we 
derive the star formation history and stretch from host 
galaxy spectrum and SN light curve, respectively. We 
then discuss the division of SNe la into a high stretch 
and a low stretch subgroup, and finally introduce the 
average spectra of Type la supernova hosts. 

3.1. A star formation history in three age bins 

We use VESPA (jToieiro et al.ll2007D to model the ob- 
served spectrum of a galaxy as a linear combination of 
up to sixteen single stellar populations (SSPs) of different 
ages and metallicities shielded by a common dust screen. 
The resulting ages, stellar masses formed, metallicities, 
and dust values for roughly 800, 000 individual galaxies 
in SDSS's Seventh Dat a Release fAbazaii an et al.|[2009l) 
have been published in jToieiro et al . (2009|). This pub- 
lic database currently holds only the SDSS Main Galaxy 
and LRG samples; the additional Stripe 82 galaxies used 
in this paper will be added soon. 

While VESPA recovers the star formation history in 
up to sixteen logarithmically spaced age bins, the num- 
ber and temporal width of these bins vary according to 
the quality of a given galaxy's spectrum. VESPA's anal- 
ysis and solutions are also model-dependent, p articularly 
for recent star formation (jToie iro et al.l[2009h . The re- 
covered mass of stars formed in the most recent age bins, 
more recently than 70 Myr ago (in the rest-frame of the 
galaxy) with our choice of bins, is particularly sensitive 
to the dust modeling. The mass formed in the next group 
of age bins, between 70 and 420 Myr, is se nsitive to the 
choic e of SSP (see Figures 15 and 20 of iTo jeiro et al.l 
l2009f ). To deal with these limitations, we: 

• select the stellar models that give the most physi- 
cally reasonable answer, and 

• limit our recovered star formation histories to three 
age bins. 

Over a typical ensemble of galaxies, the recovered star 
formation rate from 70 to 420 Myr is anticorrelated with 
and significantly lower than the rate up to 70 Myr. The 
difficulty of determining the abundance of a few hun- 
dred Myr old stars in the presence of very recent star 
formation is less pronounced with the iMarastoiil ()2005[ ) 
SSP models. To limit the impact of this, we combine all 
age bins younger than 420 Myr into a single bin. With 
these choices, the recovered star formation rate becomes 
a smooth function of lookback time. We seek the star 
formation histories of all galaxies in the same bins, and 
therefore degrade the age bins older than 420 Myr into 
two broad bins. We thus recover the star formation his- 
tories for all galaxies in the three age bins detailed in 
Table [1] For each galaxy, we obtain the stellar mass 
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TABLE \T\ 
The three VESPA age bins. 



Bin 


Age Ranj 


;e (Gyr) 


MS Spectral Types'' 


1 


0.002 - 


- 0.42 


O and B 


2 


0.42 - 


- 2.4 


A 


3 


2.4 - 


14 


F and later 



^ The main sequence (MS) spectral types that 
dominate the galaxy spectra in this age range. 



formed over each age range, the mass-weighted metal- 
licity in each bin, and an average value for the optical 
depth of interstellar dust. 

The boundaries between the age bins in Table [T] are 
chosen for convenience and to correspond to the main se- 
quence lifetimes of stellar spectral classes: we take them 
as given and estimate stellar masses. This does not imply 
that the exact age boundaries have any physical signifi- 
cance; they are simply logarithmically spaced divisions. 

3.2. Stretch measurements 

We have refitted all the light curves from SD SS-SN 
with the SALT II algorithm (|Guv et al.ll2007l ). This 
template-based routine takes the redshift and all mea- 
sured fluxes in several bands as input, and performs a 
four-parameter fit. We fit to the SDSS-SN light curves 
in the g, r, and i bands. We exclude the u band data from 
all analyses because of its poor signal-to-noise ratio, and 
exclude the z band data because Type la near-infrared 
light curves arc often double-peaked and poorly fit by 
the SALT II templates. The S ALT II outputs consist 
of a stretch- like parameter (xi in iGuv et all 2007f). colo r 
c {B — V dui B band maximum relative to Guv et all 's 
sample average), rest- frame peak B magnitude vtlb, and 
time of maximum light. For this work, we convert xi 
b ack to the t r aditio nal stretch s using the relation given 
in IGuv et al.l ((200l : 

s = 0.98 -I- 0.091x1 -I- 0.003a;? - 0.00075a;?. (3) 

Stretch s is a dimensionless parameter that describes the 
temporal width of the Type la light curve relative to a 
fiducial average. In our sample, s ranges from 0.70 to 
1.31 with a median of 0.94. The color parameter c is 
a degenerate combination of host galaxy extinction and 
intrinsic color, which for our SNe ranges from —0.14 to 
0.90 with a median of 0.06. 

While SALT II estimates the errors of all of its fitted 
parameters, we have assigned subjective quality ratings 
to each object as an independent error metric, based on 
our confidence in the best-fit value of the stretch param- 
eter s. We have excluded objects with poorly determined 
stretches (32 of the 133 with host-galaxy spectra) from 
all further analysis. Four sample light curves in g, r, 
and «, including one that was rejected in this fashion, 
are shown in Figure [T] We thus arrive at our sample of 
101 SNe selected from about 77,000 control galaxies with 
spectra. 

The selection of these 101 SNe la by inspection of the 
light curves could possibly introduce a bias in our results. 
To minimize the likelihood of this, we made all assess- 
ments of light curve quality blindly, with no knowledge 
of the host galaxy properties. Only after all of our super- 
nova light curves were fitted and rated did we perform 




Time [Days] 

Fig. 1. — Light curves in g, r, and i (top to bottom) for four 
objects with host galaxies in the SDSS spectroscopic sample; the 
bands are offset for clarity. The object shown in the upper right 
panel has little data near peak light and was therefore excluded as 
having an uncertain stretch. The other three objects entered our 
final sample of 101 SNe. Even though the object shown in the top 
left panel, 2005ex, has no data before maximum light, it is very 
well-sampled and w ell-fi t. Because we only use these fits to divide 
our sample in two ( i|3.3ll . we want to maximize our sample size by 
including SNe like 2005ex. 

any analysis on their hosts. We stress that we do not 
calibrate SNe as standard candles (which would require 
the recovery of three parameters from the light curve fit), 
but only use the light curves to divide our sample in two 
( ^3.3|) . Because of this, our light curves need not be as 
well-sampled as for cosmological analyses. We present 
details of this final sample of 101 SNe in Tabled 

3.3. The High/Low Stretch Division 

We wish to quantify the correlation between SN la 
stretch and host galaxy stellar populations using the 
VESPA analysis. Given our limited sample size, we seek 
to split our SNe into two populations based on their 
stretches. We weight each of our 101 supernovae by Mij, 
the stellar mass in its host galaxy j in a given VESPA 
age bin i. We then define the total Mi for supernovae 
with stretch s lying in the interval [a, a + As] (taking 
As — 0.04, slightly larger than a typical error in the 
stretch parameter) as 

M,(s)^ J2 M,^, (4) 

Sj G [a,a+As] 

with j running over our 101 objects. This quantifies the 
total star formation rate in a given age bin associated 
with the SNe la in our sample. Note that Mij is equal 
to the mass of possible, not actual, SN la progenitor sys- 
tems. Because galaxies are composed of a mix of stellar 
populations of different ages, we cannot derive the delay 
time distribution directly from these associated stellar 
masses. 
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TABLE [2] 
The Restricted SDSS-SN Sample 



lAU ID Host RA { J2000) Dec Stretch VESPA Host Stellar Mass (M©) 

Redshift h m s o , ii 2 - 420 Myr 0.42 - 2.4 Gyr 2.4 - 14 Gyr 



2004hz 


0, 


,1427 


02 


05 


03, 


.0 


+0 


50 


11. 


.9 


high 


8.8 X 10* 











1, 


.6 


X 


10" 


2004ia 


0, 


.1437 


02 


18 


02, 


.3 


+0 


33 


32, 


.4 


high 


7.2 X 10^'' 











2, 


.6 


X 


10" 


2004ie 


0, 


.0513 


22 


01 


46, 


.6 


+1 


14 


12, 


.0 


high 


3.7 X 10^ 


9. 


,9 


X 


10* 


2, 


.9 


X 


lOlO 


2004ig 


0, 


.1830 


00 


05 


51, 


.8 


+0 


59 


45, 


.1 


high 


3.2 X 10^ 


2. 


,9 


X 


lOlO 


1, 


.4 


X 


10" 


2005ed 


0, 


.0865 


00 


02 


49, 


.4 


+0 


45 


04, 


.3 


low 





8. 


,8 


X 


109 


4 


.0 


X 


10" 


2005ef 


0, 


.1077 


00 


58 


22, 


.9 


+0 


40 


44, 


.4 


low 


3.5 X 10** 


9. 


,5 


X 


10^ 


3, 


.4 


X 


10" 


2005eg 


0, 


.1914 


01 


02 


08, 


.5 


+0 


52 


44, 


.4 


high 


1.5 X lO'* 


4. 


,6 


X 


IQlO 


3 


.8 


X 


10" 


2005ex 


0, 


.0935 


01 


41 


51, 


.2 


+0 


52 


35, 


.0 


high 


1.8 X lO'* 











1, 


.7 


X 


10" 


2005ez 


0, 


.1298 


03 


07 


11, 


.0 


+1 


07 


10, 


.4 


low 


3.4 X 10^ 











1, 


.5 


X 


10" 


2005fa 


0, 


.1615 


01 


39 


36, 


.1 


+0 


45 


31, 


.5 


high 





1. 


,9 


X 


lOlo 


6, 


.3 


X 


10" 


2005fh 


0, 


.1190 


23 


17 


29, 


.7 


+0 


25 


45, 


.8 


high 





4. 


,9 


X 


lOlo 


2, 


.0 


X 


10" 


2005fv 


0, 


.1181 


03 


05 


22, 


.4 


+0 


51 


30, 


.1 


high 


4.4 X 10* 


4. 


,7 


X 


lOlo 


2, 


.5 


X 


10" 


2005gb 


0, 


.0864 


01 


16 


12, 


.6 


+0 


47 


31, 


.0 


low 


4.3 X 10* 


1. 


3 


X 


lOlo 


3, 


.7 


X 


10" 


2005gh 


0, 


.2591 


20 


50 


36, 


.4 


+0 


21 


14, 


.8 


high 


1.3 X 10^ 


7. 


,4 


X 


lOlo 


7, 


.5 


X 


10" 


2005gp 


0, 


.1266 


03 


41 


59, 


.3 


+0 


46 


57, 


.6 


low 


7.4 X 10* 


4. 


,6 


X 


109 


7, 


.3 


X 


lOio 


2005hc 


0, 


.0459 


01 


56 


47, 


.9 


+0 


12 


49, 


.2 


high 





7. 


,4 


X 


109 


3, 


.9 


X 


10" 


2005hj 


0, 


.0574 


01 


26 


48, 


.3 


-1 


14 


16, 


.8 


high 


9.8 X 10^ 


6. 


,0 


X 


10* 


1, 


.9 


X 


lOio 


2005hn 


0, 


.1085 


21 


57 


04, 


.2 


+0 


13 


24, 


.5 


low 


6.2 X 10* 











1, 


.1 


X 


lOio 


2005ho 


0, 


.0628 


00 


59 


24, 


.1 


+0 


00 


09, 


.4 


high 


2.3 X lO'* 











2, 


.3 


X 


lOio 


2005if 


0, 


.0670 


03 


30 


12, 


.9 


+0 


58 


28, 


.5 


low 


3.1 X 10^ 











1, 


.1 


X 


10" 


2005ij 


0, 


.1245 


03 


04 


21, 


.3 


-1 


03 


46, 


.6 


high 


2.7 X 10^ 


2. 


,0 


X 


lOlO 


1, 


.2 


X 


10" 


2005ir 


0, 


.0763 


01 


16 


43, 


.8 


+0 


47 


40, 


.4 


high 


3.1 X lO'' 


9. 


,3 


X 


109 


2, 


.9 


X 


10" 


2005je 


0, 


.0939 


02 


35 


26, 


.6 


+1 


04 


29, 


.6 


low 





1. 


,9 


X 


lOlo 


8, 


.1 


X 


10" 


2005js 


0, 


.0796 


01 


34 


41, 


.5 


+0 


36 


19, 


.3 


low 





1. 


1 


X 


lOlo 


5, 


.7 


X 


10" 


2005kt 


0, 


.0653 


01 


10 


58, 


.0 


+0 


16 


34, 


.1 


low 





1. 


,9 


X 


lOlO 


1, 


.3 


X 


10" 


2005ku 


0, 


.0454 


22 


59 


42, 


.6 


+0 


00 


49, 


.3 


high 


3.0 X lO'' 











5, 


.3 


X 


lOio 


20051j 


0, 


.0777 


01 


57 


43, 


.0 


+0 


10 


46, 


.0 


high 


5.9 X 10* 











5, 


.2 


X 


lOio 


20051k 


0, 


.1042 


21 


59 


49, 


.4 


-1 


11 


37, 


.3 


low 


7.1 X 10* 


4. 


,0 


X 


lOlo 


7, 


.9 


X 


10" 


2006eq 


0, 


.0494 


21 


28 


37, 


.1 


+1 


13 


41, 


.2 


low 


9.4 X 10^ 


8. 


3 


X 


10^ 


2, 


.1 


X 


10" 


2006er 


0, 


.0843 


00 


21 


37, 


.5 


-1 


00 


35, 


.9 


low 





1. 


,0 


X 


lOlo 


5, 


.0 


X 


10" 


2006ex 


0, 


.1472 


20 


38 


43, 


.9 


+0 


28 


28, 


.3 


high 


1.3 X 10" 


1. 


1 


X 


10" 


2, 


.7 


X 


10" 


2006fb 


0, 


.2451 


23 


35 


51, 


.5 


+0 


10 


37. 


.6 


high 


1.0 X lO'' 


4. 


,2 


X 


IQio 


3 


.6 


X 


10" 


2006fd 


0, 


.0799 


20 


37 


53, 


.2 


+1 


13 


16, 


.1 


low 





2. 


,2 


X 


lOio 


2, 


.6 


X 


10" 


2006ff 


0, 


.2353 


00 


26 


35, 


.6 


+0 


18 


07, 


.5 


high 


2.4 X lO'' 


6. 


,6 


X 


10^ 


3 


.2 


X 


10" 


2006fi 


0, 


.2306 


22 


19 


50, 


.3 


+0 


01 


27, 


.8 


high 


4.1 X lO'' 


5. 


,8 


X 


lOio 


2, 


.3 


X 


10" 


200611 


0, 


.1717 


22 


11 


27, 


.7 


+0 


45 


21, 


.5 


high 


8.7 X 10* 


9. 


1 


X 


109 


6, 


.2 


X 


10" 


2006fs 


0, 


.0991 


21 


09 


59, 


.0 


+0 


24 


31, 


.6 


high 


2.3 X lO'J 


2. 


1 


X 


lOio 


4, 


.1 


X 


10" 


2006fu 


0, 


.1985 


23 


51 


08, 


.4 


+0 


44 


46, 


.9 


high 


1.5 X lO'* 


2. 


,5 


X 


10* 


4, 


.0 


X 


lOio 


2006fv 


0, 


.1319 


01 


21 


37, 


.9 


+0 


24 


52, 


.2 


low 





8. 


,9 


X 


109 


5, 


.0 


X 


10" 


2006fy 


0, 


.0827 


23 


26 


40, 


.2 


+0 


50 


24, 


.9 


high 


2.5 X lO'* 


1. 


,2 


X 


10^ 


3 


.5 


X 


10" 


2006fz 


0, 


.1047 


00 


16 


41, 


.4 


+0 


25 


28, 


.3 


low 





2. 


1 


X 


lOio 


9, 


.6 


X 


10" 


2006gb 


0, 


.2660 


23 


59 


16, 


.5 


-1 


15 


01, 


.3 


low 


4.5 X 10^ 


6. 


,5 


X 


lOio 











2006gx 


0, 


.1807 


02 


48 


14, 


.1 


+0 


20 


49, 


.3 


high 


6.5 X lO'' 


8. 


,7 


X 


10^ 


5, 


.3 


X 


10" 


2006hd 


0, 


.2983 


21 


44 


03, 


.5 


+0 


43 


34, 


.6 


high 


2.0 X lO'' 


2. 


,9 


X 


lOio 


6, 


.7 


X 


10" 


2006hh 


0, 


.2374 


02 


42 


27, 


.0 


+0 


47 


38, 


.9 


low 


2.9 X 10* 


6. 


,9 


X 


lOio 


1, 


.2 


X 


IOI2 


2006hr 


0, 


.1576 


01 


50 


15, 


.6 


+0 


53 


14, 


.1 


high 


1.4 X 10^'' 











8, 


.4 


X 


10" 


2006hw 


0, 


.1394 


03 


13 


03, 


.4 


+0 


28 


17, 


.9 


high 


2.4 X lOi" 











3, 


.9 


X 


10" 


2006hx 


0, 


.0454 


01 


13 


57, 


.3 


+0 


22 


18, 


.0 


high 


1.3 X 10* 


1. 


,6 


X 


lOio 


1, 


.6 


X 


10" 


2006ia 


0, 


.1749 


02 


07 


19, 


.2 


+1 


15 


07, 


.5 


low 


1.5 X 10* 


4. 


,5 


X 


lOio 


1, 


.9 


X 


1012 


2006ib 


0, 


.1811 


03 


16 


11, 


.8 


+0 


36 


03, 


.4 


low 














1 


.4 


X 


10" 


2006ju 


0, 


.1486 


23 


24 


39 


.0 


+0 


43 


06, 


.0 


low 


3.5 X 10^ 


4. 


1 


X 


lOio 


8 


.4 


X 


10" 


2006jw 


0, 


.2495 


02 


23 


22, 


.3 


+0 


49 


08, 


.4 


high 


9.7 X 10* 


8. 


,1 


X 


lOio 


2 


.4 


X 


10" 


2006jz 


0, 


.1994 


00 


11 


24, 


.8 


+0 


42 


09, 


.8 


low 





1. 


1 


X 


10" 


1 


.7 


X 


1012 


2006kd 


0, 


.1363 


01 


07 


50, 


.0 


+0 


49 


41, 


.5 


high 


7.7 X 10* 











1, 


.2 


X 


10" 


2006kq 


0, 


.1983 


21 


15 


36, 


.6 


+0 


19 


17, 


.1 


high 


2.0 X 10^ 


4. 


,0 


X 


lOio 


1, 


.4 


X 


10" 


2006kw 


0, 


.1854 


02 


14 


58, 


.0 


+0 


36 


09, 


.0 


high 


9.6 X 10^ 











6, 


.2 


X 


10" 


2006kx 


0, 


.1599 


03 


42 


14, 


.7 


+0 


28 


41, 


.8 


high 


4.8 X 10^ 


8. 


,2 


X 


109 


5, 


.5 


X 


10" 


20061b 


0, 


.1819 


03 


19 


28, 


.2 


+0 


19 


04, 


.9 


high 


3.0 X 10^ 


8. 


,3 


X 


109 


5, 


.1 


X 


10" 


2006nd 


0, 


.1288 


22 


44 


59, 


.1 


-1 


00 


23, 


.8 


high 


2.4 X 10^'' 











5, 


.5 


X 


10" 


2006ne 


0, 


.0466 


01 


13 


37, 


.8 


+0 


25 


25, 


.9 


low 


1.6 X 10^ 


1. 


1 


X 


lOi" 


1, 


.7 


X 


10" 


2006ni 


0, 


.1750 


20 


54 


52, 


.4 


+0 


11 


41, 


.4 


low 





6. 


,0 


X 


lOi" 


9, 


.3 


X 


10" 


2006nn 


0, 


.1969 


01 


45 


41, 


.0 


-1 


03 


15, 


.8 


high 


1.3 X lO'* 











1, 


.3 


X 


10" 


2006nz 


0, 


.0381 


00 


56 


29, 


.2 


-1 


13 


36, 


.1 


low 





1. 


,6 


X 


109 


7, 


.6 


X 


10" 


2006oa 


0, 


.0625 


21 


23 


42, 


.9 


+0 


50 


36, 


.5 


high 





1. 


,7 


X 


10^ 


8, 


.8 


X 


107 


2006ob 


0, 


.0592 


01 


51 


48, 


.1 


+0 


15 


48, 


.3 


low 





2. 


,1 


X 


10" 


8, 


.8 


X 


10" 


2OO60I 


0, 


.1191 


23 


28 


07, 


.2 


+0 


51 


22, 


.9 


low 





3. 


3 


X 


lOio 


1, 


.1 


X 


1012 
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TABLE [2]- Continued 



lAU ID Host RA ( J2000) Dec Stretch VESPA Host Stellar Mass (Mq) 





Redshift 


h m s 




o 




// 




2 


- 


420 Myr 


0.42 - 2.4 Gyr 


2.4 - 


14 Gyr 


zUUDon 


U.U ( ly 


21 


55 


58 





— i 


U4 


1 O 7 


high 


U 








y.o X lU 


/I 7 \y 
4. / X 


iU 


zuuoop 


U.Uo41 


21 


21 


31 


o 
y 


n 


oy 


oo.y 


low 


o 
z 


A 
4 


X 


10» 


Q 7 ^ 1 n9 
o. / X lu 


O Q V 
z.o X 


lU 


zUUope 


U.ioli 


DO 


23 


09 


z 


n 




1 Q 1 

io. i 


low 


U 








o.i X iU 


i.i X 


1 ni2 
iU 


zuuDpy 


U.Uo i o 


22 


41 


42 


rt 
U 




OQ 
Uo 


1 O O 

iz.y 


low 


u 








Q 1 ^ 1 n9 

O.l X iU 


o ^ 
y.o X 


iU 


ZUU / lit 


n 0797 


00 


34 


33 


Q 
O 


1 

— i 




uo. 1 


low 


1 
i 


4 


X 


10« 


i.D X iU 


S V 
D.O X 


mil 


zuu / IIX 


n n7QA 
u.u / 


02 


06 


27 


1 
1 


-rU 




oo.o 


high 





4 


X 


lOS 


i.O X iU 


'^ A V 

0.4 X 


mil 


zuu / iiy 




03 


39 


42 


q 
o 


1 1 


Uo 


oz.z 


high 


u 








Z.Z X iU 


O.o X 


1 nil 


zuu / IIZ 


u. ioyo 


21 


03 


08 


o 

y 


— i 


Ui 


40.1 


high 


7 
i 


1 
1 


X 


109 


U 


7 V 
D. / X 


1 nil 


2007ia 


0.1310 


03 


43 


10 


1 


+0 


06 


08.9 


low 


1 


6 


X 


10^ 


O.U X iU 


5.1 X 


1 nil 
iU 


2007id 


0.1603 


21 


46 


00 


5 


-1 


13 


03.9 


high 











4.2 X 109 


1.8 X 


10" 


2007ie 


0.0934 


22 


17 


36 


7 


+0 


36 


48.0 


low 


7 


9 


X 


10« 










2007jk 


0.1829 


02 


55 


05 


6 


+0 


08 


50.8 


high 


3 


7 


X 


10* 


6.3 X 109 


1.3 X 


10" 


2007js 


0.1464 


20 


36 


48 


7 


+0 


05 


54.4 


high 


9 


4 


X 


10" 










2007jt 


0.1447 


02 


28 


32 


8 


-1 


02 


31.6 


low 


1 


7 


X 


109 


1.8 X 10^9 


9.0 X 


1019 


2007kl 


0.2571 


02 


44 


50 


9 


+0 


21 


53.4 


high 


9 


6 


X 


10* 


5.7 X 10^9 


8.5 X 


10" 


2007kv 


0.3295 


01 


10 


15 


8 


+0 


28 


19.3 


high 


3 


6 


X 


10* 


3.3 X 10^9 


3.0 X 


10" 


20071r 


0.1562 


00 


49 


00 


3 


+0 


19 


26.4 


high 


6 


1 


X 


10* 


2.9 X 10^9 


7.3 X 


10" 


2007ma 


0.1073 


00 


44 


53 


8 


+0 


59 


49.3 


high 


1 


6 


X 


109 


6.5 X 109 


5.4 X 


IQio 


2007mh 


0.1278 


03 


14 


31 


8 


+0 


16 


11.4 


high 


8 


9 


X 


10* 





4.5 X 


IQio 


2007mi 


0.1322 


03 


23 


31 


5 


+0 


39 


60.0 


low 











1.6 X 10^9 


5.2 X 


10" 


2007mj 


0.1232 


03 


34 


44 


4 


+0 


21 


19.9 


high 


2 


5 


X 


10» 


2.5 X 10^9 


1.5 X 


10" 


2007mm 


0.0665 


01 


05 


46 


7 


+0 


45 


31.8 


low 











1.5 X 109 


6.8 X 


IQi" 


2007mn 


0.0769 


02 


05 


04 





+0 


10 


28.4 


high 


2 


8 


X 


lO'J 


3.7 X 10^9 


4.4 X 


10" 


2007nj 


0.1540 


02 


52 


27 


4 


+0 


15 


06.6 


low 


4 


7 


X 


109 


1.7 X 10^9 


3.8 X 


10" 


2007ok 


0.1655 


02 


28 


24 


3 


+0 


11 


04.8 


high 


3 





X 


10* 


3.6 X 10^9 


4.3 X 


10" 


2007ol 


0.0560 


01 


37 


23 


7 


+0 


18 


43.2 


low 











7.4 X 109 


1.7 X 


10" 


2007om 


0.1052 


23 


54 


20 


7 


+0 


55 


03.4 


high 


3 


9 


X 


10» 


7.9 X 10^9 


7.6 X 


10" 


2007OU 


0.1132 


02 


23 


42 


7 


+0 


49 


33.6 


high 


2 


2 


X 


109 





1.2 X 


10" 


2007ph 


0.1294 


20 


51 


13 


4 


+0 


57 


20.9 


low 











1.1 X 10^9 


9.2 X 


10" 


2007pt 


0.1752 


02 


07 


38 


5 


+0 


19 


26.4 


high 


1 


8 


X 


109 





7.6 X 


10"' 


2007px 


0.1080 


00 


22 


44 





+0 


28 


44.4 


high 


2 


6 


X 


1019 





5.3 X 


10" 


2007py 


0.2094 


03 


29 


31 


6 


+0 


30 


56.0 


low 


1 


6 


X 


109 





1.3 X 


10" 


2007qa 


0.1085 


01 


52 


33 


9 


+1 


14 


38.7 


high 


1 


6 


X 


109 





1.3 X 


10" 


2007qr 


0.1359 


02 


52 


29 


2 


-1 


08 


22.3 


high 


6 





X 


10» 


1.4 X 109 


5.7 X 


10"' 


2007rs 


0.1241 


00 


46 


27 


4 


-1 


03 


44.1 


low 











6.6 X 10^9 


1.8 X 


1012 



Position of SN; data available from_ http://www.cfa. harvard. edu/iau/lists/Supernovae.html^ and 
|http://www.cfa.harvard.edu/iau/cbat.htrQl| 
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With this caveat in mind, the resulting distributions 
of stellar mass as a function of stretch in different mass 
bins are shown in Figure [2] The top panel is a histogram 
of stretches, while the second, third, and fourth panels 
weight each SN by its host stellar mass in young, middle, 
and old age bins, respectively. The second and fourth 
panels, in particular, show that lower stretch SNe live 
in hosts with large quantities of old stars but little re- 
cent star formation, while higher stretch SNe are in hosts 
with abundant recent star formation. Given the rough 
bimodality shown here, we call Type la supernovae high 
stretch if they have s > 0.92, and low stretch otherwise. 
In SDSS-SN, this division yields 60 high stretch and 41 
low stretch SNe; the binnings are Usted in Table [2l 

This resu l t conf irms the w ork oflSullivan et al.l ( 2006[ ) , 
INeiU et all (pOOOt ). GaUaghereLalJ (12005D, al^ others: 
that high stretch SNe are associated with star-forming 
galaxies. We now seek to study this association in more 
detail and to derive a quantitative delay time distribution 
for high and low stretch SNe la. 

3.4. Average Spectra 

The difference in stellar populations indicated by Fig- 
ure [2] should manifest itself in the spectra of the hosts. 
To examine this, we constructed average spectra of the 
high stretch and low stretch hosts. We corrected the 
spectra for Mi l ky Wa y exti nction using the dust maps of 
ISchlegel et all (|1998f ) and a iFitzpatricH (|T999D Ry = 3.1 
extinction law. We scaled each spectrum to a common 
r band fiber magnitude, shifted the spectra to z = 0, 
masked bad pixels flagged by the SDSS pipeline, and 
weighted by the inverse variance in coadding. To obtain 
robust estimates of the mean and variance of the aver- 
age spectra, we have used bootstrap resampling on the 
samples of 60 high stretch and 41 low stretch hosts. 

The results, shown in Figure [31 are striking. The aver- 
age spectrum of a high stretch host shows exceptionally 
strong nebular emission lines such as the Balmer series, 
[Oiii], [On], [Nil], and [Sii] and a strong blue contin- 
uum, all clearly indicative of recent and ongoing star 
formation. The average spectrum of low stretch super- 
nova hosts, by contrast, shows absorption lines charac- 
teristic of old stellar populations and little evidence of 
interstellar gas. The difference spectrum (lower panel) 
looks much like the spectrum of a main-sequence B type 
star with nebular emission lines superimposed. Intrigu- 
ingly, a single B3 spectrum from Silva & Cornell (1992) 
fits the continuum and absorption lines of the difference 
spectrum better than models of either a 50 Myr-old single 
stellar population, or 400 Myr of continuous star forma- 
tion. While the bootstrap errors are too large to exclude 
either of these possibilities. Figure [3] tantalizingly sug- 
gests very young progenitors for some SNe la, as the dif- 
ference spectrum appears to require a population < 50 
Myr old (B3 stars live for about 30 Myr). Note that 
we cannot be absolutely certain that the progenitors are 
this young, if bursts of star formation typically last long 
enough so that a somewhat older population is likely to 
be associated with ongoing star formation. 

Figure [3] confirms the difference in host galaxies indi- 
cated in Figure [5] by the VESPA-derived star formation 
histories. High stretch, luminous Type la's are associ- 
ated with young O and B stars, while lower stretch SNe 
are found in galaxies with much older stellar populations. 




Stretch s 



Fig. 2. — Histograms of stellar mass of supernovae la hosts binned 
by SN stretch. Objects in the top panel are unweighted, while SNe 
in the lower three panels are weighted by the stellar mass of their 
host galaxy in young (top), middle-aged (middle) and old (bottom) 
stars. SNe with s > 0.92 (to the right of the dotted line) are 
associated with young stars, while lower stretch SNe are not. To 
further explore this distinction, we define supernovae with s > 0.92 
as high stretch SNe la. 

Figure [3] offers a clear and dramatic confirmation of the 
association between high stretch SNe and star forma- 
tion. The average spectra change imperceptibly when 
the high/low stretch boundary is varied from 0.90 to 0.92, 
representing a typical stretch error as reported by SALT 
II. Our results are therefore unlikely to be affected by 
measurement errors in the stretch parameter. 

4. TOWARDS A DTD 

While Figure[3]is compelling on its own, we can do bet- 
ter. With a well-controlled survey, spectra for all galaxies 
in the control sample, and VESPA star formation histo- 
ries for all of these galaxies, we now seek to calculate a 
delay time distribution in the three age bins of Table [TJ 
We do not assume any functional form for the efficiency 
of making SNe as a function of progenitor age, but simply 
see what constraints the data alone can provide. 

Our method is as follows: 
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High Stretch Hosts 
Low Stretch Hosts 




Bootstrap Mean Spectra 



Bootstrap la Errors 

rMM L 



JLjl. 



— Smoothed Difference 
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4000 5000 6000 7000 
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8000 



Fig. 3. — Coadded SDSS spectra of hosts of 60 high stretch 
(s > 0.92 - blue) and 41 low stretch {s < 0.92 - red) supernovae 
(upper panel); the shading shows the la error range derived from 
bootstrap resampling. Note the strong blue continuum and nebular 
emission lines, indicative of recent and ongoing star formation, in 
the high stretch hosts. The low stretch host spectrum has features 
characteristic of much older stellar populations. The difference 
(lower panel) looks remarkably similar to the spectrum of a typical 
B star (note the strong Balmer series in absorption), suggesting 
young progenitors for high stretch Type la SNe. The stellar spec- 
trum is from Silva & Cornell (1992). 



1. Parametrize the DTD as an explosion rate per unit 
stellar mass for stars in each of three age bins i, 

2. Assume a prior probability distribution on the 
DTD explosion rates, 

3. Select DTDs from this prior, 

4. Generate samples of mock hosts from the DTDs 
and VESPA star formation histories of our control 
sample, 

5. Compare the average spectra of the mock hosts and 
observed hosts to calculate a likelihood for each 
DTD realization, 

6. Repeat steps 3-5 many times to explore the Ci 
space and obtain a posterior probability distribu- 
tion on the DTD. 

We now discuss the general form of the delay time dis- 
tribution and then each of these steps. 



4.1. The Delay Time Distribution 

Our star formation histories consist of the total stel- 
lar mass formed in each of three age bins for each host 
galaxy. The most general delay time distribution is there- 
fore a set of three efficiencies, e^, representing the mean 
number of SNe per unit stellar mass per year from pro- 
genitors in age bin i (Table [1]) . We treat high stretch 
and low stretch SNe separately in this analysis, giving 
two sets of efhciencies eh,i and 6;.^. 

We may effectively remove two of these six parameters 
by requiring that the sets of efficiencies be appropriately 
normalized, i.e. that the total number of expected SNe 
times their probability of detection equal the number of 
SNe observed, 



NsNM = ^ ph,j{detect)t^^ 

gals j 



(5) 



where phj{detect) is the probability that a supernova of 
high stretch in galaxy j enters the sample, Mij is the 
amount of stellar mass formed in age bin i in galaxy j, 
NsNji is the number of high stretch SNe la observed, and 
t is the duration of the survey. An identical constraint 
applies to low stretch supernovae. 

For future use, we also introduce normalized explosion 
efficiencies corresponding to the fraction of high or low 
stretch supernovae produced by stars in a given age bin, 
e.g. 



1 



N. 



SNJi 



Ph,j (detect) Mi^ 



(6) 



gals j 



Thus, ejj j and ej ^ each sum to unity. 

4.2. Priors 

In order to use a Bayesian analysis to constrain the €h,i 
an d e;.i, we need t o choose a prior distribution. Follow- 
ing |L^laci ()1814[ ) and general practice in the literature, 
we seek a uniform prior. However, due to the normaliza- 
tion constraint, we cannot place uniform priors simulta- 
neously on all of the eh,i or e;_i. An alternative would be 
to choose our priors to be uniform in area over the two 
dimensional plane defined by Equation ([5]). However, the 
projection of this prior onto any of the e'^ j yields a prob- 
ability distribution strongly peaked towards low values. 
Physically, this biases us against a DTD in which one 
age bin is responsible for a large fraction of the SNe (see 
Figure H]). To avoid these combinatorial effects, we use 
a Monte Carlo sampling that places a uniform prior on 
one of the e'^ ^ and one of the e[ ^ in each realization. 
By choosing i randomly in each iteration, we retain the 
symmetry between all of the ej^ ^ and ej ^ . 

We note that there are many possible priors that look 
qualitatively similar to that shown in Figure |4l We do 
not argue that our choice is optimal, only that it is rea- 
sonable. We have also tested the dependence of our con- 
straints on the choice of prior and found little variation, 
as long the prior probability density does not approach 
zero for any values of e[. 

4.3. Creating a Mock Sample 

Given a model delay time distribution (i.e., sets of eh,i 
and ei^i), we compute the expected number of high and 
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Fig. 4. — Projected prior distributions for each e^, the fraction of 
high/low stretch SNe with progenitors in age bin i. The normaUzed 
e'^ are related to the physical explosion rates by Equation A 
prior uniform in the 2-D plane defined by Equation (O would be 
strongly biased against DTDs concentrated in one age bin. 

low stretch supernovae (<S 1) observed in each galaxy in 
the control sample. For a galaxy j, this is the sum of the 
explosion efficiencies times the mass in each stellar age 
bin weighted by the detection probability: 



{nsN,h,]) = Phjidetect) ^ eh,iMi^j. 



(7) 



Our code thus generates realizations of {e'^} for both high 
and low stretch, converts to {ei} using Equation (jH]), and 
computes the number of SNe la in each control galaxy as 
a Poisson random number with a mean given by Equation 

©■ 

For the detection probability, we use a crude estimate 
of the selection function and tweak the parameters to 
match the approximate redshift range of the subset of 
SDSS-SN with host galaxy spectra. We adopt the func- 
tional form 



p{detect) 



(8) 



where rripeak is the supernova's peak SDSS r band mag- 
nitude, miim is an approximate limiting magnitude of 
the survey, and a is a softening parameter to account for 
the dispersion in SN peak magnitudes and colors and the 
survey detection efficiency. Because this model does not 
naturally capture the fact that higher stretch SNe are 
above any threshold for a longer period of time, we al- 
low the parameters and a to differ for high and low 
stretch objects. Fitting to the observed redshift distribu- 
tion, we adopt — 19.5 and a — 0.6 for low stretch 
SNe, and nium = 20.2 and a = 0.4 for high stretch SNe. 
The peak supernova magnitude in r is given by 

mpeak = Mstand ~ a{s - 1) + fl{z) + K[s, z) + Ar (9) 

where M stand is the standard ized peak B absolute magni- 
tude, a is the iPhillipsI ()1993D parameter relating stretch 
and peak luminosity, s is a typical stretch, is the dis- 
tance modulus, Ar is the r band Milky Way extinction 
along the supernova's line of sight fiSchlegel et al.,il998i ). 



Fig. 5. — Stretch vs. color {B — V at peak B magnitude relative to 
a fiducial average), as determined by SALT II. The two quantities 
are not correlated, which allows us to fold the color variation into 
the detection efficiency (Equation [Sj. 

and K{s, z) is an app roximate JT-co rrection to the SDSS 
r band at 2; = (see iNugent et al.l (j2002) for details on 
the JT-correcti on). We adopt th e values M^tand — —19.41 
and a = 1.56 ()Guv et al.ll2005[ ). and take the typical val- 
ues of s to be our sample medians, 1.02 for high stretch 
and 0.83 for low stretch SNe la. Since we fit for 
and a to reproduce the observed redshift distribution af- 
ter fixing all other parameters, changing the standard- 
ized absolute magnitude Mp^ak will not affect our re- 
sults. Because we fit high and low stretch SNe sepa- 
rately, the value of a similarly has no effect. Equation 
([9]) does not account for variation in color (i? — ^ at peak 
B magnitude), which is instead absorbed into the soft- 
ening parameter a. If color is correlated with stretch. 
Equation ([9]) will not fully capture differences between 
the detectability of high and low stretch SNe. However, 
as shown in Figure [SJ we find no evidence of such a cor- 
relation. 

Figure [B] compares the redshift distribution of a mock 
sample generated by Equations ([7]), ([8]), and ([9]) with 
the observed redshift distribution. While there are slight 
differences, the means and widths of the distributions 
agree reasonably well for both high and low stretch SNe. 
The detection probability given by Equation (jS]) is in- 
tended only to be approximate and to set the redshift 
range of galaxies that serve as potential hosts. Modest 
variations in Equations ([8]) and ^ have almost no effect 
on our results, particularly because we compare spectra 
weighted by their inverse variances: the slight discrep- 
ancy at high redshift is dominated by galaxies with poor 
signal-to-noise ratio spectra that have a very small im- 
pact on the average spectrum. As an additional check, 
we have verified that Equations (O, ([8]), and (HJ success- 
fully reproduce the observed range of host galaxy masses 
as derived from VESPA. 

4.4. A Likelihood Function 

To use Bayes' Theorem, we need to calculate the likeli- 
hood of observed host properties given a sample of mock 
hosts selected with a specific model for the {ci}. We use 
a comparison metric based on the average spectra of the 
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Fig. 6. — Redshift histograms of the observed SNe with SDSS 
host galaxy spectra against the histograms of one realization of 
mock hosts with the detection probability given by Equation 
We adopt values of mjim = 20.2 and a = 0.4 for high stretch SNe, 
^Um = 19-5 and a = 0.6 for low stretch SNe. 

host galaxies; essentially, with our Monte Carlo realiza- 
tions of delay time distributions, we seek to reproduce 
the mean spectra shown in Figure [3l We use a statis- 
tic to compare the average spectra of our Monte Carlo 
realizations to those of the observed hosts. 

To obtain useful constraints, we need both very high 
signal-to-noise ratio average spectra for the mock sam- 
ples, and reliable estimates of the mean and variance as- 
sociated with the average spectra of the observed hosts. 
The first requirement is met by using a large mock sam- 
ple drawn from our 77,000 control galaxies. The second 
problem is solved using bootstrap resampling of the 60 
high stretch and 41 low stretch actual hosts, as described 
in §3.41 and shown in Figure [3l As a check, we compared 
the bootstrap errors to the error-weighted pixel-by-pixel 
scatter between spectra and found very good agreement. 

With an average spectrum and bootstrap errors for 
the SN hosts, we may assign a goodness-of-fit value 
to the mean spectrum of the mock hosts by averaging 
the values over all pixels. We then use this reduced 
to compute a rejection probability. 

There are two reasons why we have not directly com- 
pared the stellar populations of SN la hosts and non- 
hosts. First, such a comparison would require us to ap- 
ply VESPA to the SN hosts to derive their star formation 
histories. We are limited by ^/N statistics in the hosts, 
and errors in the VESPA outputs are much larger than in 
the spectra. Second, the recovered stellar masses would 
be used twice for the control galaxies: once to determine 
how many supernovae to set off in the mock samples 
(Equation [7]) , and once to compare the resulting mock 
hosts to the observed hosts. Any metric would therefore 
contain the sum of the squares of the recovered stellar 
masses, so that random variances in these masses would 
add and produce a bias. Our chosen metric, by combin- 
ing spectra rather than derived quantities, avoids both 
drawbacks. 



5. RESULTS 

We have run 100,000 Monte Carlo realizations of the 
delay time distribution from the prior shown in Figure 21 
In each realization, we draw a set of normalized explo- 
sion efficiencies (e'. Equation [5]) , convert to the physical 
explosion efficiencies (e. Equation [5]) , and use Equations 
(O, ([5]), and ^ to generate a sample of mock hosts. We 
then construct the composite spectra for the mock high 
stretch and low stretch hosts and compare them to the 
spectra of observed hosts shown in Figure [3] using a 
test. We take the likelihood of each realization to be the 
rejection probability (Q function) of its computed value 
of x^- In practice, our best- fit models do produce for- 
mally good fits, with x^ per pixel typically in the range 
0.65 - 1.05. Normalizing these likelihoods by their sum 
over all Monte Carlo realizations, we obtain the poste- 
rior probability distributions for the explosion efficien- 
cies. The distribution for a single efficiency is computed 
by integrating over the other variables. 

The resulting probability distributions are shown in 
Figure [Tj The lower horizontal axes give the explosion 
efficiencies in physical units {e.h,i and eu, Equation [5]) , 
while the upper axes show the fraction of all SNe la in the 
SDSS spectroscopic sample formed as high or low stretch 
from a given progenitor age bin. These are computed 
from the eh,i or e/^i as, e.g., 

fh,i = 4r^ Ph,jidetect)M^j, (10) 

gals J 

where ph,j{detect) is the detection probability for high 
stretch SNe in galaxy j, Mi^j is the stellar mass formed 
in age bin i in galaxy j, NgN is the total (high and low 
stretch) number of SNe detected, and the sum is taken 
over control galaxies. Equation (llOp differs from Equa- 
tion ([6]) only in the denominator: while ej^ ^ and ej ^ each 
sum to unity, fh.i and fi^i together sum to unity. Note 
that because a volume-limited sample of galaxies will 
have different properties from the SDSS spectroscopic 
sample (in particular lower masses), fh.i or for the 
full sample of SDSS supernovae may differ from the val- 
ues we recover. 

These same quantities, / and e, are also plotted in Fig- 
ure [51 with the top panel showing the supernova rate in 
physical units {th,i in blue and e/,i in red) and the bottom 
panel showing the fraction of SNe from a given range of 
progenitor ages {fh,i in blue and fi^i in red). The dots 
show the median values, while the colored and black bars 
represent the 68% and 95% confidence intervals, respec- 
tively. For distributions peaked near zero, only upper 
limits and one-sided confidence intervals are shown. The 
constraints on / and e are also listed in Tabled 

The results are striking. We can constrain most high 
stretch SNe la to have progenitors younger than < 400 
Myr, and Figure[3]tantalizingly suggests an even younger 
characteristic age (B3 stars live for ~30 Myr). While 
young stars dominate the production of high stretch 
SNe la, they make no significant contribution to the low 
stretch rate. Instead, these supernovae have a character- 
istic delay time of at least 2-3 Gyr. Intermediate pro- 
genitors, evolving on the time scale of a Main Sequence 
A star, contribute little to either the high or low stretch 
rate. In the SDSS Stripe 82 galaxies, the prompt, high 
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TABLE [3] 

The DTD: results of the Monte Carlo 



Fraction of Total Type la SNR f 
0.1 0.2 0.3 



Age Ranj 


;e (Gyr) 


eh 


f b 
Jh 




f b 

Jt 


0.002 - 
0.42 - 

2.4 - 


- 0.42 

- 2.4 
14 


'0-16 

< 13 
< 0.48 


'^■^0-0.05 

< 0.20 

< 0.18 


< 13 

< 2.2 

1 cro + O.Oe 

i.as_Q 12 


< 0.04 

< 0.04 

c:o+0.02 
'J-OS-O.OS 



Note. — The errors given are 95% confidence intervals. 
Defined in Equation (0, units are 10"^* SNeMg^yr'^. 
^ Defined in Equation JTOj. 

stretch channel and the delayed, low stretch channel each 
account for roughly half of all Type la supernovae. 

The formal constraints we obtain for the low stretch 
rate from young stars (Figure [7]) are especially stringent. 
Because even blue galaxies have lots of old stars, we find 
it difficult to reproduce a spectrum as quiescent as that 
shown in Figure [S] a sample selected purely by old stellar 
mass will still include star-forming galaxies. Our mea- 
sured spectra are only consistent with a delay time dis- 
tribution in which all (or nearly all) low stretch SNe la 
are produced by old stars. Nevertheless, we cannot rule 
out any progenitor channel expected to produce < 1 de- 
tectable supernova over the entire population of SDSS 
galaxies, and we have therefore set a floor on our SN 
rate constraints from Poisson statistics. This floor de- 
termines the limits on e;_i and e;^2, the low stretch rates 
from young and middle-aged stars, in Table[3]and Figure 

m 

To compare our results with earlier work based on 
the A + B model (Equation [T]), we need to interpret 
A and B in our framework. We take A, the coeffi- 
cient of the total stellar mass, to be our supernova rate 
for old stars, £3. We convert B, the coefficient of cur- 
rent star formation (in SNe Mq^), to our rate for young 
stars €1 (in SNeMg^yr"^) by multiplying by the tem- 
poral width of our first age bin, about 400 Myr. Our 
result s agree reasonably well with those of iNeill et aLl 
(poM ) (young ei = 200 ± 50, old eg = 1.4 ± 1.0 in our 
u nits of 10"^"^ SNeMq Vr-^), but less well with those 
of ISullivan et~aI1 (12006^ ( ei = 100 ± 20, 63 = 5.3 ± 1.0) 
or iMannucci et"all (|2005h (d = 650 ± 280 or SOOti^^, 

ea = 4.4ti;4). Much o f the discrepan cy i s due to the defi- 
nition of stellar mass: ISullivan et al.l and IMannucci et al.l 
apply a correction for dead stars, d ecreasing stell a r mas s 
and correspondingly increasing e. IDildav et al.l ([2008h . 
who measured these parameters using SDSS-SN with- 
out correcting for dead stars, obtained ei = 230l;7g ^^^"^ 
£3 — 2.8 ± 1.2 . In addition, the values of depend on 
the authors' choice of stellar models and proxies for re- 
cent star formation, and may therefore be expected to 
disagree. Further, while our relative rates for progeni- 
tors of different ages are extremely robust to the details 
of the detection function (Equation[S]), the normalization 
of our rates is less so. A poor estimate of the detection 
function would multiply all of the rates eh,i and ei^i in 
Table [3] by a constant of order unity. 

6. DISCUSSION AND CONCLUSIONS 
6.1. Robustness and Systematics 
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Fig. 7. — The posterior probability distributions of the DTD. The 
lower horizontal axes give the explosion efficiencies ^ and e; i, 
the supernova rates per unit stellar mass per year for progenitors 
in age bin i (Equation [Sj. The upper axes show j and /; j, 
the proportions of the curr ent observed SN la rate in the SDSS 
Stripe 82 galaxies (Equation [TOj. Nearly all high stretch SNe have 
progenitors < 400 Myr old (top panel), while Figure [3l suggests 
that the typical age may be even younger. Low stretch SNe have a 
characteristic delay time of at least 2-3 Gyr (bottom panel), with 
essentially no contribution from young progenitors. 
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Fig. 8. — The recovered DTDs, expressed in physical units (Equa- 
tion \5\ to p p anel) and as a proportion of the total SN la rate 
(Equation 1101 bottom panel) in the SDSS spectroscopic sample. 
Upper limits are shown for distributions peaking near zero, medi- 
ans and two-sided confidence intervals are shown otherwise. Points 
are drawn at the midpoint of their relevant age ranges (see Table 
[T]|, and the high and low stretch bars are offset for ease of viewing. 

SDSS-SN is a large, controlled, untargeted survey, and 
is therefore largely free of systematics in target selection. 
The subset with host galaxy spectra does have additional 
selection criteria that can lead to systematics; however, 
the host galaxies were selected from the entire Stripe 
82 spectroscopic sample only through the occurrence of 
a detectable SN la. Further, our comparison of Monte 
Carlo and observed host galaxies is based not on derived 
quantities (like stellar masses), but rather on the spec- 
tra themselves. We therefore believe our results to be 
robust. Here, we address several possible sources of bias 
and argue that all should be minor. 

One possible source of systematics is VESPA and its in- 
put stellar models (see H3.ip . We have chosen our models 
and the resolution of our recovered star formation histo- 
ries to minimize these effects. More importantly, we have 
avoided any use of VESPA on the actual hosts (other 
than in Figure [U which we use solely to set the stage) , 
and only use it to select our samples of mock hosts. 
Therefore, the star formation histories from VESPA will 



only bias our results if they are incorrect in the mean. 
The clear evidence of young stars in Figure [3] qualita- 
tively supports the results in FigurejTl 

Our earlier paper, lAubourg et all ([20Q8l) . used VESPA 
to show the existence of a prom pt component as sociated 
with the age bin to 180 Myr. lAubourg et al.l and the 
present manuscript use different approaches to recover 
information about SN la progenitors. While in this paper 
we aim at recovering the full DTD in the least paramet- 
ric possible way allowed by the data, in lAubourg et al.l 
we focused on recovering the shortest age of the prompt 
progenitor channel. It is possible that very young star 
formation can mask the presence of older stars. This 
in itself does not remove the requirement for the pres- 
ence of a young population, as implied by the analysis of 
lAubourg et al.l but in order to avoid this issue, we have 
taken a cautious approach and combined the young bins 
into a broad bin. However, the difference spectrum in 
Figure [3] suggests that the youngest bin in our study may 
be dominated by progenitors substantiall y younger than 
400 Myr, which would give support to the lAubourg et all 
findings. Note that we cannot with these data exclude 
the possibility that the progenitors are older, as we have 
demonstrated only that very recent star formation is re- 
quired to be present. An older progenitor population 
could be reconciled provided correlations in the star for- 
mation rate mean that older star formation is, in the 
majority of cases, accompanied by a very young popula- 
tion. We are limited by the signal-to-noise ratio of the 
difference spectrum; with a larger supernova host sam- 
ple, we should be able to model better the age of the 
stellar populations of the difference spectrum. 

Systematic errors in the shape of the observed spec- 
tra will give errors in the VESPA-derived star formation 
histories. In a sample of physically identical galaxies, 
such an effect could, for example, produce a subsample 
with bluer measured spectra and thus younger inferred 
stellar populations. Selecting galaxies based on their de- 
rived masses of young stars will therefore produce aver- 
age spectra that are too blue. In this way, differences 
in the average spectra of mock hosts selected by Equa- 
tion d?]) will be exaggerated by using the VESPA-derived 
masses rather than the (unknown) physical masses. How- 
ever, this effect would make it easier to reproduce the 
differences in average spectra seen in Figure [3l biasing 
us against the different progenitor ages for high and low 
stretch SNe la seen in Figures [71 and El These effects are 
also small: lAdelman-McCarthv et alT (|2008[ ) show that 
spectrophotometric calibrations are accurate to about 
4% rms (see their Figures 4 and 5). 

Another possible source of systematics is the variation 
in SN and host galaxy properties. Galaxies with young 
stellar populations tend to have more dust, which could 
render SNe associated with younger stellar populations 
fainter. We have not taken host galaxy dust into account 
other than by adding scatter to our selection function, 
through the a parameter in Equation ([8]). However, the 
colors of SNe are measurable and depend in part on host 
galaxy extinction. As shown in Figure El the color pa- 
rameter recovered by SALT II is uncorrelated with the 
stretch. Perhaps more surprisingly, it is also uncorre- 
lated with the mean interstellar dust extinction as fit 
by VESPA. Any systematic variations of color with host 
galaxy properties must therefore be small, and should 
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h ave little effect on our results. 

iKellv et al.l ()2009D have found a hint of another sys- 
tematic effect, a correlation between peak SN la lumi- 
nosity and host galaxy mass not captured by the color or 
stretch variation. However, the magnitude of this effect 
is too small (< 0.1 magnitude) to significantly impact 
our study. 

In addition to these systematics, we could suffer bi- 
ases either from our light curve ratings fi )3.2p or from 
our choice of parameters in Equations ([5]) and To 
eliminate the former, we did all of the light curve fits 
and ratings blindly (! j3.2p . with no knowledge of the host 
galaxies. For the l atter, we took aU free parameters in 
Equation ^ from lGuv et~all (|2005D and fit the param- 
eters miim and a in our detection efficiency (Equation 
[5]) to match the observed redshift distributions of both 
high and low stretch SNe. Changing and a could 

multiply all of the recovered explosion efficiencies e by 
a constant of order unity, but our results are otherwise 
insensitive to the details of the selection function. 

Finally, we measure the SN la rates in the SDSS 
spectroscopic sample, which contains fewer low-mass, 
metal-poor galaxies than a volume-limited sample would. 
Should SN la pr oduction be suppressed in low metallicity 
environments as iKobavashi et al.l (|1998f) and others have 
suggested, SN la rates could be lower in a volume- limited 
sample than our results would indicate. 

6.2. Connection to Progenitor Models 

It is remarkable how cleanly our observational cut 
in stretch divides the sample into groups with distinct 
progenitors. The choice of 0.92 as a division is some- 
wha t arbitrary: it was a guess motivated by Figure [2] 
and ISullivan et ahl ()2006l ). The observed distribution of 
stretches does not show a clear bimodality (see Figure 
[S]) and indeed appears continuous. It is of course possi- 
ble that SNe la have more than two progenitor channels, 
each distinguished by a range of stretches. However, Fig- 
ure [2] provides little guidance on where to split either the 
high or low stretch sample, and in any case, our sample 
sizes of 60 high stretch and 41 low stretch SNe are too 
small to profitably subdivide further. 

Given the division by stretch into two samples, we can 
test the compatibility of our DTD with predictions of var- 
ious progenitor models. We find that many can explain 
the prompt channel but have extreme difficulty repro- 
ducing our low stretch DTD, in which nearly all systems 
take more than ^2.4 Gyr to explod e. The theoretical s in- 
gle degenerate DTD calculated by iWang et all (|2010[ ) is 
compatible with our high stretch sample but falls short 
of our low stretch rates by at least an order of magni- 
tude (see their Figures 8 and 9). AU of the DTDs pub- 
lished by Greggio' (2005) have too high a rate from young 
stars to match our low stretch sample, th ough several 
are co mpatible with our high stretch rates. iRuiter et al.l 
()2009() calculate two-peaked theoretical DTDs for vari- 
ous common envelope scenarios, but with the exception 
of one semidetached double white dwarf binary model, 
find rates in young stars inconsistent with our low stretch 
constraints. It is possible that a single progenitor chan- 
nel could produce luminous, high stretch objects with 
a short delay time and subluminous, low stretch objects 
with a long delay time, perhaps because of different com- 
positions at white dwarf birth. Though recent progress 



has been made (jWooslev et al.|[2007l ). we still lack an un- 
derstanding of the dependence of SN la stretch on the 
progenitor composition, and do not attempt to address 
the likelihood of this scenario in this paper. 

6.3. Conclusions 

We have used SDSS-SN, the only large untargeted su- 
pernova survey at low redshift, to constrain the progen- 
itor populations of Type la supernovae. The blind na- 
ture of the survey renders it free of most systematics, 
while our use of spectra of both hosts and a large control 
sample has allowed us to minimize the impact of stellar 
population models on our results. We dramati cally con- 
firm the two populations seen bv ISullivan et a l. (2006), 
iMannucci et al.' (2006D, iScannapieco fc Bildsteni ()2005i) . 
Aubourg ct al. (200^ and others, finding a "prompt" 
component of luminous, high stretch SNe with a char- 
acteristic delay time < 400 Myr (the difference spectrum 
in Figure [3] hints at a time as short as tens of Myr) and 
a "delayed" component of subluminous, low stretch SNe 
with a delay time > 2.4 Gyr. While our results are in 
broad agreement with the A + B model (Equation [T]) , 
they place strong constraints on the progenitor ages and 
cause difficulties for many extant theoretical delay time 
distributions. We caution against any physical interpre- 
tation of the precise values of the age boundaries, which 
serve only as limits of integration. 

Type la supernovae have played a key role in our 
current understanding of the cosmological model. In 
spite of our incomplete understanding of Type la pro- 
genitors and explo sions, SN la s urveys like SNLS 
(lAstier et al.l l2006l) . ESSENCE (IW;ood-Vasev et al.l 
120071) and SDSS-SN (IKessler et al.l 120091) continue 
to provide the best constraints on cosmological pa- 
rameters. Further improvements, for example with 
the JDEM candidate SuperNova Acceleratio n Probe 
(|SNAP Collaboration: G. Aldering et alll200l . wih re- 
quire systematic effects in SN la standardization to be 
controlled to 1 — 2%. The identification of the stellar 
evolution paths that can yield SNe la is therefore a key 
issue: the brightness of a supernova could depend on the 
nature of its progenitor, and the demographics of SNe 
la would depend on redshift through the evolution time. 
This could yield an effective, non-cosmological, depen- 
dence of the mean "standardized" absolute magnitude on 
redshift that would bias dark energy measurements. In 
principle such an effect could be advantageous, as analy- 
ses of host spectra could be used to determine the prob- 
abilities of each progenitor rou te, and hence the use of 
the appropriate iPhillipsI (|1993D stretch-luminosity rela- 
tion. This will require both better constraints on Type 
la progenitor models and an improved understanding of 
the Phillips relation. 

While Type la supernovae are most widely studied be- 
cause of their use as standard candles, they are also dy- 
namically important in the interstellar medium and are 
believed to be the main source of iron-peak elements. 
A fast route to Type la supernovae therefore has im- 
plications for the interpretation of alpha-enhancement. 
If most SNe la are old, iron enrichment will be signif- 
icantly delayed from the onset of star formation, and 
alpha-enhancement will be associated with short-lived 
periods of relatively recent star formation. Significant 
early iron production by "prompt" SNe la would mod- 
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ify this conclusion, making it more difficult to explain 
alpha- enhancement and supporting alternative explana- 
tions such as a modification of the initial mass function. 

Future studies will certainly improve on our results; 
the Baryon Oscillation Spectroscopic Survey (BOSS, 
ISchlegel et b1\\2007 ). part of SDSS-III^, is perhaps the 
most promising. BOSS is taking spectra of all super- 
nova hosts from SDSS-SN, expanding our sample by a 
factor of ~3. While a control sample is not naturally de- 
fined for these SN hosts, the untargeted nature of SDSS- 
SN means that a control group could be built out of 
a deep, volume-lim ited sample of galaxies like GAMA 
(jBaldrv et al.l 120091 ). Such a sample could be used to 
significantly tighten our error bars. 

Increasing the temporal resolution of our recovered 
DTD may be more challenging. This will require better 
spectra, both of SN hosts and control galaxies, and es- 
pecially better stellar population synthesis models. The 
latter point is crucial: we had to combine all ages younger 
than 420 Myr into a single bin because of biases intro- 
duced by the modeling of dust and stellar spectra. Split- 
ting this age bin will allow us to better determine the 
timescale of "prompt" SNe la and thus the initial masses 
of their progenitors. 

With better stellar population synthesis models, 
BOSS, GAMA, and other surveys would hold tremen- 
dous promise. They would allow us to construct a sam- 
ple several times as large as the one used here, and with 
better recovered star formation histories. With such a 
sample, we could improve both the precision and tem- 
poral resolution of our DTD by a factor of ^-^2, placing 
strict constraints on theoretical SN la progenitor models. 
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